Drude weight and total optical weight in a t-t'-J model 
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We study the Drude weight D and the total optical weight K for a t-t'-J model on a square lattice 
that exhibits a metallic phase-modulated antiferromagnetic ground state close to half-filling. Within 
a suitable 1/N expansion that includes leading quantum-fluctuation effects, D and K are found to 
increase linearly with small hole doping away from the Mott metal-insulator transition point at 
half-filling. The slow zero-sound velocity near the latter transition identifies with the velocity of 
the lower-energy branch of the twofold excitation spectrum. At higher doping values, D and K 
eventually saturate and then start to decrease. These features are in qualitative agreement with 
optical conductivity measurements in doped antiferromagnets. 
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I. INTRODUCTION 

The optical properties of the metallic state that 
emerges upon doping the insulating Heisenberg antifer- 
romagnet with mobile holes are of great interest because 
such a system, currently described by a t-t'-J model, is 
believed to capture low-energy physics of the high-T c su- 
perconducting copper-oxide layers Q. In the context of 
the latter model we have recently calculated using 
a suitable 1/N expansion, the real part of the finite- 
frequency optical conductivity, i.e., the optical absorp- 
tion cr(u>), and found that its magnitude scales with small 
hole concentration while its low-frequency line shape dis- 
plays a structureless broadband regime. These results 
are in qualitative agreement with the generic properties 
of the midinfrared band observed in optical conductiv- 
ity measurements |^-|^] in doped antiferromagnets. Here 
we extend our earlier work by presenting results for 
the Drude weight D and the total optical weight K that 
determine the zero- frequency response of the system and 
the optical conductivity sum rule, respectively. As it is 
well-known Q , the Drude weight D serves as an order pa- 
rameter for the Mott metal-insulator transition occurring 
in such a system, hence the significance of its dependence 
on hole concentration and relevant coupling constants. 
Furthermore, the fraction of the total optical weight re- 
siding at zero frequency, i.e., the fraction D/K, defines 
the inverse of the mass enhancement factor. The latter 
quantity is a direct measure of the strength of the quasi- 
particle interactions or equivalently, in the context of the 
1 /N expansion, of the strength of the quantum fluctua- 
tions. 

The t-t'-J model Hamiltonian can be expressed in 
terms of the Hubbard operators \ ah = \a)(b\ as 
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spin-down electron, and the summation convention is in- 
voked. Here J is the antiferromagnetic spin-exchange in- 
teraction between nearest-neighbor sites on a square 
lattice, while for the hopping matrix elements Uj we as- 
sume 



t if i,j are nearest neighbors, 
tij = { —t' if i,j are next nearest neighbors, 
otherwise. 



(2) 



Our conventions in (|^) incorporate opposite signs for t 
and t' as dictated by quantum-chemistry calculations 
|t],H for Cu-0 clusters and fits of the shape of the Fermi 
surface observed by angle-resolved photoemission spec- 
troscopy ||. In Ref. |lCj we generalized the local con- 
straint associated with (111) to \i° + xf = N, where N 
is an arbitrary integer, and simplified the commutation 
properties of the x ab 's to be those of the generators of the 
U(3) algebra. A generalized Holstein-Primakoff realiza- 
tion of this algebra in terms of Bose operators can then 
be employed to develop a perturbation theory based on 
the 1/N expansion, restoring the relevant physical value 
N = 1 at the end of the calculation. For an average 
electronic density n e close to half-filling (n e < 1) and a 
sufficiently large t', the ground state of (|l|) in the large- 
N limit is an unconventional metallic antiferromagnetic 
(AF) state characterized by the usual (tt, tt) modulation 
of the spin configuration and an unusual (tt, — tt) phase 
modulation of the condensate |l(| . It is the optical prop- 
erties of this phase-modulated AF state that the present 
work is concerned with. 

The remainder of this paper is organized as follows. 
Section [0] contains a brief summary of the relevant Kubo 
formalism for the optical conductivity. In Sec. Ill we de- 



rive the analytic expressions of the 1/N expansion for 
the Drude weight and the total optical weight, includ- 
ing leading quantum-fluctuation effects. 
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where the index corresponds to a hole, the Greek indices 
/X, v, . . . assume two distinct values, for a spin-up and a 



In Sec. 

present explicit numerical results for the latter quantities 
along with some further comments on the optical absorp- 
tion and a discussion of the zero-sound velocity which is 
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related to the Drude weight and the inverse compressibil- 
ity. Our conclusions are summarized in Sec. 



II. KUBO FORMALISM 

Using standard linear-response theory, the Kubo for- 
mula for the real part of the frequency-dependent optical 
conductivity cr(cj) at zero temperature is expressed as the 
sum of two physically distinct terms P~0|l2| 



Re[a(cu)} = Tre 2 D5(uj) + ct(lu) 



(3) 



The first term in (^|) , involving the delta function at zero 
frequency (li = 0), is due to the free acceleration of the 
charge-carrying mass by the electric field and therefore 
the associated spectral weight, i.e., the Drude weight D, 
should vanish in the insulating state and be finite in the 
metallic state ||. D measures the ratio of the density 
of the mobile charge carriers to their mass. The second 
term in (||), called the optical absorption cr(w), is due 
to finite-frequency {lj > 0) optical transitions to excited 
quasiparticle states. More explicitly we have that 



= -t E I H J l°) \ 2s ^ - E o) <"] 
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where the summation is taken over a complete set of en- 
ergy eigenstates |m) with excitation energies (E m — Eo) 
above the ground state |0), A is the total number of lat- 
tice sites, and J is one of the Cartesian components of 
the current operator J associated with (hi), 
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(5) 



where R 4 is the position vector for site i. Similarly, the 
Drude weight D is given by 



n=K 2 v |(m|J|0)| 2 
e 2 A 2-, E m -E ' 

where K is the expectation value 

K=j£(0\-T\0) , 
with the operator T defined as 

T = -^t ij |R i -R,| 2 X f X f , 



(6) 



(7) 



(8) 



z — 4 being the coordination number of the square lat- 
tice. Integrating now over frequency both sides of (^), 
using the identity f °° dui S(lo) = | and the explicit forms 
(§ and (§), we arrive at the well-known optical conduc- 
tivity sum rule p3[ 



jf dwRe[B{u)] = —K , 



(9) 



which, indeed, identifies K as the total optical spectral 
weight. 

It should be noted that in lattice models involving only 
nearest-neighbor hopping, T defined in (0) coincides with 
the kinetic energy operator given by the first term on the 
right-hand-side of (Q). However, the presence of an ad- 
ditional next-nearest-neighbor hopping, i.e., a hopping t' 
along the diagonal of the square unit cell as in the t-t'-J 
model under study, invalidates the latter simple result 
noting that: |Rj — Rj| 2 = 2, for next-nearest-neighbor 
sites Therefore in this more general case, the total 
optical weight K given by (|?]) is not just the kinetic- 
energy expectation value in the ground state. Never- 
theless, by an elementary application of the Hellmann- 
Feynman theorem K can still be extracted from the 
ground-state energy (0Lff|0) by differentiation with re- 
spect to the hoppings [H t and t': 



(10) 



The identity (|l7j) is actually implemented in the following 
section to derive the leading terms of the 1 /N expansion 
of K from the corresponding terms of (0|i? |0) which our 
theory readily provides. 



III. 1/N EXPANSION 

The 1/N expansion of the Hamiltonian (|l|) and the 
optical absorption (Q) around the phase-modulated AF 
configuration, up to and including terms of order N, has 
been carried out in our earlier works |h]J^]- For ease of 
reference, we quote in ([TT|)— (p~4|) the main results neces- 
sary for our present development. Specifically, the ex- 
pansion of the Hamiltonian (TP) reads pi] 



H = N 2 AE a + NAEi 



(11) 



where Eg and E\ are the classical (large-iV limit) and 
zero-point energy per lattice site, respectively, for the 
relevant physical value N = 1: 

zJ 

E a = -zt'n e (l - n e ) —n 2 , 

(12) 

Ex = -^-( 2 ~ ne ) ~ + 2A~ E t Wl ( q ) + W2 ( q )] ' 



In the above Aq and -B q are the normal-mode operators 
that diagonalize the leading (quadratic) quantum fluctu- 
ations while w„(q), n = 1,2, are the dispersions of the 
two branches of the spectrum of elementary excitations. 
An immediate consequence of (|ll]) is the following result 
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for the l/N expansion of the ground-state energy, up to 
and including terms of order N, 



|#|0) = N 2 AE + NAEx . 



(13) 



The corresponding expansion of the optical absorption 
M reads U 



cr(w) = N y ' x 



ujA 



q 



2 (q),5[cj 1 (q)+tJ 2 (q)-w] . (14) 



The analytic expressions for the dimensionless function 
w 2 (q), corresponding to the current matrix elements, and 
the dispersions uj n (q), n = 1,2, are rather involved and 
are summarized in the Appendix. The result (Jl4|) has 
been extensively discussed in our earlier work of Ref. |^] 
while some further comments are added in Se c. |TV| . 

At present, an important point to note in (|l4|) is the 
lack of a classical contribution, i.e., a term of order A^ 2 , 
to <t(lu). In other words, the finite- frequency component 
<j(lu) of the optical conductivity is due exclusively to the 
quantum fluctuations whose leading contribution, within 
the present l/N expansion scheme, is of order N. This 
conclusion is consistent with similar observations made 
Bang and Kotliar jL5| for the finite- frequency optical con- 
ductivity of the simple t-J model treated by a slave-boson 
diagrammatic l/N expansion technique. 

The derivation of the l/N expansion of the Drudc 
weight D and the total optical weight K, up to and in- 
cluding terms of order N, proceeds now in two steps: 

(i) The result (13) is used to implement identity (p^[). 

(ii) The second term on the right-hand-side of (||) is ex- 
pressed initially as the integral of over positive fre- 
quencies and then ( |l4| ) is exploited to obtain immediately 
its l/N expansion. Taking into account (|l2|), the final 
analytic results of this straightforward procedure may be 
written in the form 



D = N 2 D Q + ND 1 , 
K = N 2 K + NK! , 



(15) 
(16) 



where Dq and K$ are the classical (large- AT limit) contri- 
butions to the Drude weight and the total optical weight, 
respectively, 



D = K 



2/5 , 5 \ 



(17) 



while D\ and K± are the corresponding contributions due 
to leading quantum- fluctuation effects, 



A wi(q) + w 2 (q) 



(18) 
(19) 



The equality D = K quoted in (|l7|) is a consequence 
of the absence of quantum fluctuations in the large-A^ 
limit and implies that in this classical approximation, de- 
scribed by the terms of order A^ 2 , the 5-function Drude 
peak centered at zero frequency (u = 0) carries the total 
optical weight. However, as noted earlier on, the lead- 
ing quantum fluctuations, described by the terms of or- 
der A", already give rise to optical absorption and hence 
distribute part of the total optical weight to finite fre- 
quencies (u> > 0). The presence of the terms of order 
A" in ([l5]) and ( |l6| ) leads therefore to D < K; a generic 
property anticipated also from (^|) and (|l9|). 

In view of ( |l2| ) and the explicit forms quoted in the Ap- 
pendix, the derivatives with respect to t and t' appearing 
in (|l^) can be carried out analytically. The remaining 
wave vector integrations over the square Brillouin zone 
in ( |l8| ) and (|l9|), as well as in jl^), can then be per- 
formed numerically Explicit numerical results derived 
from (|l5|)-(|l9|) in the way just described are presented in 
the following section. 



IV. EXPLICIT RESULTS 

For completeness and better insight to the results for 
the Drude weight and the total optical weight that will 
be presented shortly, it is instructive to begin our discus- 
sion here with some examples of the optical absorption 
line shape a(u>). In Fig. [l] we draw cf{lo), determined 
from (|l4|) with N = 1, for typical values of the dimen- 
sionless ratios e — t'/t, t/J and the hole concentration 
(l — n e ) that are relevant for the copper-oxide layers 
As noted in Ref. the limiting value a(u — > 0) is fi- 
nite, at finite hole doping, while for small doping val- 
ues (such as the 10% hole doping considered in Fig. [I]) 
the high-frequency divergence of <r(u>) comes from the 
q = (it/2, 7t/2) point of the Brillouin zone. The frequency 
position of this peak (divergence) is then determined as 

z J 

[wi(q) + w 2 (q)] q =( 7r / 2 ,7r/2) = zt'(2 - n e ) + —n e , (20) 

and apparently is independent of the nearest-neighbor 
hopping t. For the parameter values corresponding to the 
solid line of Fig. |](a), the peak position of <t(lu) calcu- 
lated from (^0|) is given by u> = 3. 78 J. The latter value, 
with an estimated Jl6|-[l8| J « 0.13-0.16 eV, accounts 
aptly for the 10 w 0.5 eV peak of the midinfrared band 
observed in optical conductivity measurements |^-^| in 
doped antiferromagnets. 

From the explicit form ( p0| ) follows that close to half- 
filling and for parameters satisfying Its/ J < 1 the peak of 
<j(oj) shifts slowly to lower frequencies upon hole doping, 
i.e., with increasing (1 — n e ). Such a systematic shift of 
the peak of the midinfrared band is, indeed, observed in 
the afore- mentioned optical conductivity measurements. 
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Interesting enough, close to half- filling the same condi- 
tion 2te/J < I has been shown in Ref. [jlO| to ensure 
the softening of the magnon velocity upon hole doping; a 
trend that is also observed in inelastic neutron-scattering 
measurements from doped antiferromagnets |l6| , |r^ ] . The 
present theory provides therefore a context for the quali- 
tative understanding of both experimental observations. 
Concluding our comments on cr(ui) we note that (|3|) and 
(|) yield J™ducr(L)) = ?f(K-D). Therefore, the 
spectral weight carried by the finite-frequency compo- 
nent a(u>) is given by the difference (K — D) and hence 
its dependence on the parameters e, tj J, and n e follows 
from that of D and K. 

Let us now consider the large-TV limit contributions 
Dq and Kq to the Drude weight and the total optical 
weight, respectively. The equality of these two quan- 
tities has already been discussed following (19). Here 



we note that the vanishing overlap between the opposite 
sublattice spin states in the AF configuration, along with 
the absence of quantum fluctuations in the large- N limit, 
leaves the direct hopping t' between same sublattice sites 
as the only relevant process of charge transport in this 
classical approximation. This argument makes plausible 
the independence of D , K from t and J seen in (|l7|). 

Furthermore, the n e -dependence in ( |l7| ) implies that: 
(i) Do = Kq oc n e , for n e — > 0, and (ii) Do — Kq oc 
(1 — n e ), for n e — > 1. The linear increase of Do and 
Kq with small electron density n e away from the empty- 
lattice limit (n e = 0) is, of course, an expected behavior 
in this "free" electron gas regime. On the other hand, the 
linear increase of Do and Kq with small hole concentra- 
tion (1 — n e ) away from the half-fillcd-band limit (n e = 1) 
is an important consequence of the local constraint which 
prohibits the occupancy of any lattice site by more than 
one electron and leads inevitably to the Mott insulator 
at half-filling (D = = Kq, for n e = 1.) The latter 
property, typified in point-(ii) above, serves to interpret 
the "free" charge carriers near half-filling as being holes 
rather than electrons. Then, the position of the maxi- 
mum of Do located at quarter-filling (n e — i) provides 
an estimate of the critical amount of doping at which the 
character of the charge carriers changes from holelike to 
electronlike with increasing (1 — n e ). The present large- 
TV limit results capture already generic features of the 
n e -dependence of the Drude weight and the total opti- 
cal weight that are derived by exact diagonalization of 
the simple t-J model and the large-?/ Hubbard model on 
small clusters |f^,^,|l] . The reason can be traced to the 
fact that our generalized Holstein-Primakoff realization 
p0[ resolves explicitly the local constraint and thus the 
ensuing 1 /N expansion incorporates already part of the 
strong correlation effects, implied by this constraint, in 
the leading order. 

As a check of consistency of the quantum-liquid de- 
scription emerging from the present 1/JV expansion 



scheme it is worthwhile at this point to discuss briefly 
the zero-sound velocity, i.e., the velocity of the long- 
wavelength charge excitations, in the physically relevant 
regime near half- filling (n e — > I). As demonstrated in de- 
tail in Ref. [[t0| , a separation of spin and charge degrees 
of freedom sets in progressively and becomes asymptot- 
ically exact as the half-fillcd-band limit is approached. 
The latter limit permits then an accurate classification 
of the elementary excitations in terms of a mode that 
describes the spin excitations and a mode that describes 
the charge excitations. At long-wavelengths (|q| — > 0), 
in particular, the latter mode corresponds to the lower- 
energy branch of the twofold excitation spectrum given 
by W2(q) = c 2 |q|. Thus near half- filling [n e — > I) 
the velocity of the long-wavelength charge excitations is 
given by c 2 . On the other hand, in the context of con- 
ventional quantum-liquid theory p^j the zero-sound ve- 
locity may be expressed in terms of the Drude weight 
Do and the inverse compressibility 1/k = tv^Eq , where 
E' Q ' = d 2 E Q /dn 2 e , as y/D /(n*K) = yj ' D E' ' . Hence the 
consistency of the quantum-liquid description requires 
the validity of the identity: c 2 = UDqEq , for n e — > 1. 
The asymptotic form of c 2 for n e — > 1 has been derived 
in Ref. [[t0| by considering explicitly the long-wavelength 
limit of w 2 (q). Here, using ( |l2] ) and ( |l7| ) to calculate 
yj DqEq and comparing with the latter result for c 2 one 
can easily verify that indeed p3] 



c 2 = VDoE'J = 2e 



for 



"I 1/2 



1. (21) 



In view of the radically different way in which c 2 and 



•J DqEq are derived, the relation (|2lj) provides a strin- 
gent consistency check of our theory and serves to iden- 
tify unambiguously c 2 as the slow zero-sound velocity for 
n e — > 1. Furthermore, from ( |l2| ) follows that Eq is in- 
dependent of n e . Hence (|2l| ) reveals that the physical 
reason for the softening of the velocity c 2 of the lower- 
energy branch w 2 near the Mott metal-insulator transi- 
tion at half-filling (c 2 oc y/l — n e — * 0, for n e — * 1) is 
the collapse of the Drude weight which in turn is a direct 
consequence of the local constraint. It should be noted 
that this slow zero-sound mode near the Mott transi- 
tion appears also in slave-boson treatments of the local 
constraint when the fluctuations of the associated sta- 
tistical gauge field are taken into account beyond the 
mean- field approximation [ p3[ . Concluding our discus- 
sion of the zero-sound velocity we note that away from 
the asymptotic regime n e —> 1 and with increasing hole 
concentration, w 2 (as well as uj\) involves an increasingly 
strong hybridization between spin and charge degrees of 
freedom and consequently c 2 corresponds no longer to 
the velocity of the long- wavelength charge excitations. 
Having analyzed the classical (large- TV limit) contribu- 
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tions to the Drude weight and the total optical weight, 
we consider in the remainder of this section the complete 
results for D and K that include the leading quantum- 
fluctuation corrections, as determined from ( |l5| ) and (16), 
respectively, with N = 1. In Fig. ^ (for e = 0.45) and 
Fig. U (for e = 0.40) we draw as a function of the hole 
concentration (1 — n e ): (a) the Drude weight D, (b) the 
total optical weight K, and (c) the fraction D/K, for 
t/J = 1.0 (solid line) or t/J = 1.5 (dashed line). For 
comparison, in each figure we also depict by dotted line 
the classical (large- N limit) result for the drawn quan- 
tity: (a) Do, (b) Kq, and (c) Dq/Kq. We remind that 
the latter classical results, determined from (|l7|), are in- 
dependent of the ratio t/J. Our main observations from 
Fig. H and Fig. || are summarized as follows. 

(i) D and K increase linearly with small hole concen- 
tration (1 — n e ) away from the Mott metal-insulator tran- 
sition point at half- filling (n e = 1), consistent with the 
notion that the "free" charge carriers in this regime are 
the holes. However, as a result of the leading quantum 
fluctuations, D and K are no longer equal. Close to half- 
filling and for the typical values of e and t/J considered, 
the fraction D/K is reduced from its "classical" value of 
1 to about 0.5 and is almost independent of (1 — n e ) in the 
doping range of interest; see the solid and dashed lines 
in Fig. H|(c) and Fig. ||(c). These results are consistent 
with optical conductivity measurements J^-|{| in doped 
antifcrromagnets. In particular, for e = 0.45, t/J = 1.0, 
corresponding to the solid line of Fig. 0(c) , and a typical 
10% hole doping we have: D/K = 0.52, corresponding 
to a mass enhancement factor K/D = 1.92. The latter 
value compares reasonably well with the experimental 
estimate |] (K/D) cxp w 2.3. Furthermore, our typical 
value D/K = 0.52 is compatible with the corresponding 
prediction of 0.6 derived by exact diagonalization |2Q| 
and anyon techniques p4]j for the simple t-J model, and 
implies that 52% of the total optical weight resides at 
the zero-frequency Drude peak ire 2 D5(u>) while the rest 
is carried by the finite-frequency component cr(u>). It 
should be noted that the strength of the quasiparticle 
interactions, as inferred from the mass enhancement fac- 
tor K/D rss 2.0, is not exceptionally large B. This fact, 
however, is not surprising noting that important correla- 
tions effects, implied by the local constraint, are already 
absorbed into the charge carrier density which near half- 
filling is given by the hole concentration (1 — n e ) that 
vanishes at n e = 1, leading to the Mott metal-insulator 
transition. 

(ii) At higher doping values, D and K eventually sat- 
urate and then start to decrease. However, the leading 
quantum fluctuations shift the maximum of D and K 
away from its "classical" location at quarter-filling to- 
wards half-filling. This shift is small and slightly differ- 
ent for D and K. As noted earlier on and emphasized 
by Dagotto et al. [ pp[ , the position of the maximum of 
D provides an estimate of the doping value at which the 



charge carriers change from holelike to electronlike with 
increasing (1— n e ), hence at roughly the same position the 
Hall coefficient Rh should change sign from positive to 
negative. For e = 0.45 (0.40) and t/J = 1.0, correspond- 
ing to the solid line of Fig. |2|(a) [Fig. |^(a)], the maximum 
of D occurs at the hole concentration (1 — n e ) = 0.44 
(0.36). The latter result is then compatible with numeri- 
cal and analytical calculations |25| of Rh, in the context 
of the simple t-J model, that predict a sign change at 
about 0.40-0.33 hole concentration, and with Hall effect 
measurements [ p6| in doped antifcrromagnets that report 
a sign change of Rh at about 0.3 hole concentration. 

(iii) We cannot approach the "free" electron gas regime 
close to the empty-lattice limit (n e = 0) because the 
phase-modulated AF configuration, around which the 
present l/N expansion is carried out, becomes unstable 
beyond a critical value of the hole concentration (1 — n e ). 
The intervening instability is reflected in the elementary 
excitations whereby the velocity ci of the lower-energy 
branch u>2 becomes zero at the critical doping value jlO) . 
This results in a very rapid decrease of D, K, and D/K in 
the immediate neighborhood of the latter doping value; 
a behavior clearly seen in Fig. | and Fig. § In this 
critical doping regime, the magnitude of the (negative) 
leading quantum-fluctuation corrections D\ and K\ be- 
comes even larger than that of the classical (large-JV 
limit) terms Dq and Kq, and therefore the l/N expan- 
sion breaks down. 

(iv) Finally, by comparing the solid and dashed lines 
m Fig. | and Fig. | we conclude that at any given dop- 
ing value D, K, and D/K all increase with increasing t' 
and/or J. This trend is consistent with similar results 
derived by exact diagonalization of the t-J model, ex- 
tended to include the so-called three-site term jl9|], and 
of the Hubbard model ^o|. Therefore, this trend should 
be regarded as being generic for models involving a direct 
or an effective hopping between same sublattice sites. 



V. CONCLUSIONS 

In this paper we have presented a detailed study of the 
Drude weight D, the total optical weight K, and the frac- 
tion D/K that defines the inverse of the mass enhance- 
ment factor, in the phase-modulated AF state of the t- 
t'-J model (0)-(|). The generic features of these quanti- 
ties observed in optical conductivity measurements |^-[^] 
in doped antifcrromagnets are qualitatively reproduced 
when the leading quantum- fluctuation effects, around 
the afore-mentioned semiclassical ground state, are taken 
into account within a suitable l/N expansion. 

Specifically, D and K increase linearly with small hole 
doping away from the Mott metal-insulator transition 
at half-filling, consistent with the notion that the "free" 
charge carriers in this regime are the holes. Our theoret- 
ical prediction of a mass enhancement factor K/D w 2.0 
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that is almost independent of the hole concentration, in 
the doping range of interest, compares reasonably well 
with the corresponding experimental estimate. Q With 
increasing hole doping D (and K) reaches eventually a 
maximum at a value that is compatible with measure- 
ments [^6| of the critical doping at which the Hall coeffi- 
cient Rh changes sign from positive to negative, signaling 
a change in the character of the charge carriers from hole- 
like to electronlike. Furthermore, for typical parameter 
values the peak position ( |20| ) of the finite-frequency com- 
ponent <t(uj) of the optical conductivity accounts aptly 
for the experimentally observed 0.5 eV peak of the mid- 
infrared band d ||] . 

Finally, the slow zero-sound velocity near the Mott 
metal-insulator transition point at half-filling is shown to 
identify with the velocity of the lower-energy branch of 
the twofold excitation spectrum, thus providing a strin- 
gent consistency check of the quantum-liquid description 
emerging from the present 1/N expansion scheme. 



where for an arbitrary wave vector q 
plicit forms of the coefficients w^' 

zt' n 



qj Aq, 



(q x ,q y ) the ex- 
and Va read 



A q = 



4 (l-n e ) 
zt' M 



zJ 



(1 — COS q x cos q y ) 
■cos,q x cosq y ) , 



1 + - {cosq x + cosq y ) 



,(+) 



,(-) 



,(+) 



zJ 

— n e (cos q x 



iq y ) 



zt 



n e (cosq x 



''Qy) 



= i/ ( + ) 



+ A q - r q + zt'(l - n e ){l 



(A4) 
(A5) 
(A6) 
(A7) 

cos<j x cosq y ) , 
(A8) 
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APPENDIX: 

In this Appendix we summarize the analytic expres- 
sions for the dispersions w n (q), n = 1,2, of the two 
branches of the spectrum of elementary excitations, and 
the dimensionless function w 2 (q), corresponding to the 
current matrix elements, that enter the optical absorp- 
tion (FLJ) and the leading quantum-fluctuation correc- 
tions fl8| ) , ([l9]) for the Drude weight and the total optical 
weight. 

Following the notation conventions of Ref. (l(J, the 
dispersions w n (q), n — 1,2, are given by 



,?(q) = R(q) + 2^S(q) 



{(q) = R(q) - 2 A /5(q) 



with 
R(d) = 

and 
S(q) 



,(+) 



,(-) 



,(+) 



(Al) 



(A2) 
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^q 'q U q A q 



U q T q 



(A3) 



The analytic expression for the function it 2 (q) is very 
involved. However, the result can be presented in a com- 
pact manner with the help of a four-component matrix 
notation . To this end we introduce the 4x4 matrices 
pi, P2, and P3, defined as pi = ai <S> 1, i = 1,2,3, with 
the ai denoting the ordinary 2x2 Pauli matrices. The 
same symbols will also be used in the following to de- 
note the 4x4 matrices 1 ® cr,, i = 1, 2, 3, for notational 
simplicity. We then have 



w 2 (q) = |[«»(q)+«v(q)] 



(A10) 



where u 2 (q), a — x,y, are positive-definite dimension- 
less functions expressed as traces over 4x4 matrices. 
Specifically, 



<(q) = - 7 



1 



A J a;i(q)w 2 (q)5(q) 
+Tr[v a (q)A(q, ui (q)K (q)A(q, -w 2 (q))] 
+Tr[w a (q)A(q,w 2 (q))?; a (q)A(q,-a; 1 (q))]} , 

(All) 

where v Q (q), a = a;, y, and A(q, ui) are the 4x4 real 
matrices defined below: 



v x (q) = smq x 



( 1 - c osq y -^ ( I 



--n e cosgj,cri + -n e pi<r 2 



3n,, 



173 
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u y (q) = sin q y 



-ell 



3n e \ 1 
— I cos + - I 1 



--n e cosq x ai - - / n e p\a 2 
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(A13) 
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with 



W(q) = 4+V 3 
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and 

M{o L ) = 



C3 



q q 



^q T q 



Pi 0-2 



(A16) 



In conclusion it is worth emphasizing that the disper- 
sions w n (q), n = 1, 2, defined in (Al) as well as the func- 
tion w 2 (q) defined in ( A10 ), depend on q only through 
the factors cos q x and cos q y and are symmetric under the 
interchange of variables (q x ,q y ): Lu n (q x ,q y ) = u n (q y , q x ), 
n = 1, 2, and u 2 (q x ,q y ) = u 2 (q y , q x ). 



[22] Note that an expression analogous to yjD E'^ for the 
velocity of the long-wavelength charge excitations in the 
one- dimensional Hubbard model appears implicitly also 
in the work of H. J. Schutz, Phys. Rev. Lett. 64, 2831 
(1990). 

[23] J. P. Rodriguez, Phys. Rev. Lett. 73, 1675 (1994); Phys. 

Rev. B 44, 9582 (1991); 45, 5119(E) (1992). 
[24] A. M. Tikofsky, R. B. Laughlin, and Z. Zou, Phys. Rev. 

Lett. 69, 3670 (1992). 
[25] B. S. Shastry, B. I. Shraiman, and R. R. P. Singh, Phys. 

Rev. Lett. 70, 2004 (1993); M. Deeg, H. Fehske, and H. 

Biittner, Europhys. Lett. 26, 109 (1994); E. Dagotto, A. 

Nazarenko, and M. Boninsegni, Phys. Rev. Lett. 73, 728 

(1994). 

[26] H. Takagi et al, Phys. Rev. B 40, 2254 (1989); H. Y. 
Hwang et al, Phys. Rev. Lett. 72, 2636 (1994). 



[1] For a recent review, see E. Dagotto, Rev. Mod. Phys. 66, 
763 (1994). 

[2] G. C. Psaltakis, Phys. Rev. B 51, 2979 (1995). 

[3] S. L. Cooper et al, Phys. Rev. B 41, 11605 (1990). 

[4] J. Orenstein et al, Phys. Rev. B 42, 6342 (1990). 

[5] S. Uchida et al, Phys. Rev. B 43, 7942 (1991). 

[6] W. Kohn, Phys. Rev. 133, A171 (1964). 

[7] M. S. Hybertsen, E. B. Stechel, M. Schluter, and D. R. 
Jennison, Phys. Rev. B 41, 11068 (1990). 

[8] G. A. Sawatzky, in Earlier and Recent Aspects of Super- 
conductivity, edited by J. G. Bednorz and K. A. Miiller 
(Springer- Verlag, Berlin, 1990), p. 345. 

[9] J. Yu and A. J. Freeman, J. Phys. Chem. Solids 52, 1351 

(1991) . 

[10] G. C. Psaltakis and N. Papanicolaou, Phys. Rev. B 48, 
456 (1993). 

[11] B. S. Shastry and B. Sutherland, Phys. Rev. Lett. 65, 
243 (1990); B. S. Shastry, Mod. Phys. Lett. B 6, 1427 

(1992) . 

[12] D. J. Scalapino, S. R. White, and S. C. Zhang, Phys. Rev. 

Lett. 68, 2830 (1992); Phys. Rev. B 47, 7995 (1993). 
[13] R. A. Bari, D. Adler, and R. V. Lange, Phys. Rev. B 

2, 2898 (1970); P. F. Maldague, ibid. 16, 2437 (1977); 

D. Baeriswyl, C. Gros, and T. M. Rice, ibid. 35, 8391 

(1987). 

[14] D. Baeriswyl, J. Carmelo, and A. Luther, Phys. Rev. B 
33, 7247 (1986); 34, 8976(E) (1986). 

[15] Y. Bang and G. Kotliar, Phys. Rev. B 48, 9898 (1993). 

[16] G. Aeppli et al, Phys. Rev. Lett. 62, 2052 (1989). 

[17] J. Rossat-Mignod et al, in Dynamics of Magnetic Fluc- 
tuations in High- Temperature Superconductors, edited by 
G. Reiter, P. Horsch, and G. C. Psaltakis (Plenum, New 
York, 1991), p. 35. 

[18] S. Shamoto et al, Phys. Rev. B 48, 13817 (1993). 



FIG. 1. Optical absorption <t(lj) (in units of e 2 /K) over 
hole concentration (1 — n e ) vs frequency for 1 — n e = 0.10, 
t/J — 1.0 (solid line) or t/J = 1.5 (dashed line), and (a) 
£ = 0.45, (b) e = 0.40. 



FIG. 2. For e = 0.45 and t/J = 1.0 (solid line) or t/J = 1.5 
(dashed line), including leading quantum-fluctuation effects: 
(a) Drude weight D vs hole concentration, (b) Total optical 
weight K vs hole concentration, (c) Fraction D/K vs hole 
concentration. In each figure, the dotted line is the classical 
(large- N limit) result for the drawn quantity and is indepen- 
dent of t/J according to Eq. (p^). 

FIG. 3. For e = 0.40 and t/J = 1.0 (solid line) or t/J = 1.5 
(dashed line), including leading quantum-fluctuation effects: 
(a) Drude weight D vs hole concentration, (b) Total optical 
weight K vs hole concentration, (c) Fraction D/K vs hole 
concentration. In each figure, the dotted line is the classical 
(large- iV limit) result for the drawn quantity and is indepen- 
dent of t/J according to Eq. (|ri|). 
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